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Weak antilocalization (WAL) and linear magnetoresistance (LMR) are two most commonly observed 
magnetoresistance (MR) phenomena in topological insulators (TIs) and often attributed to the Dirac 
topological surface states (TSS). However, ambiguities exist because these phenomena could also come from 
bulk states (often carrying significant conduction in many TIs) and are observable even in non-TI materials. 
Here, we demonstrate back-gated ambipolar TI field-effect transistors in (Bi 0 .o4Sb 0 .96)2Te 3 thin films grown 
by molecular beam epitaxy on SrTi0 3 (l 11), exhibiting a large carrier density tunability (by nearly 2 orders of 
magnitude) and a metal-insulator transition in the bulk (allowing switching off the bulk conduction). 
Tuning the Fermi level from bulk band to TSS strongly enhances both the WAL (increasing the number of 
quantum coherent channels from one to peak around two) and LMR (increasing its slope by up to 10 times). 
The SS-enhanced LMR is accompanied by a strongly nonlinear Hall effect, suggesting important roles of 
charge inhomogeneity (and a related classical LMR), although existing models of LMR cannot capture all 
aspects of our data. Our systematic gate and temperature dependent magnetotransport studies provide 
deeper insights into the nature of both MR phenomena and reveal differences between bulk and TSS 
transport in TI related materials. 

Topological insulators (TIs) are an exotic state of quantum matter with nominally-insulating bulk and spin- 
momentum-locked Dirac fermion conducting surface states, promising potential applications in both 
nanoelectronics and spintronics 1 \ Many interesting phenomena based on such topological surface states 
(TSS) have been proposed, such as Majorana fermions 4,5 , exciton condensation 6 , topological magnetoelectric 
effect 7 , etc. Bi 2 Se 3 , Bi 2 Te 3 , and Sb 2 Te 3 have been proposed and identified as prototype 3D TIs, possessing TSS with 
a linear dispersion of energy vs momentum in the bulk band gap (BBG) 811 However, most of these commonly 
studied TI materials inevitably have impurities or defects during their growth, resulting in a doped bulk which can 
account for a significant part of the conductance, making it difficult to study and utilize the novel electronic 
transport of the TSS. To access the surface transport properties of the 3D TIs, various strategies have been 
attempted to suppress the bulk conduction, for example by (compensation) doping, increasing the surface-to- 
volume ratio, or electric gating. A ternary compound (Bi 1 „ x Sb x ) 2 Te 3 , an alloyed mixture of Bi 2 Te 3 and Sb 2 Te 3 , has 
been shown as a promising group of 3D TIs with excellent tunability of the electronic properties by varying the 
composition x 1213 . Recently, the quantum anomalous Hall effect 14 has been successfully observed in Cr doped 
(Bi x Sb!. x ) 2 Te 3 , adding further interests to examine the underlying electronic transport properties of this ternary 
TI system. 

Most transport studies of TIs have focused on magnetotransport 15 , particularly magnetoresistance (MR). In a 
few experiments, Shubnikov-de Haas oscillations 111618 with n Berry phase were observed and used as direct 
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transport evidence for TSS Dirac fermions. However, the majority of 
TIs do not show such quantum oscillations because of disorder or 
inhomogeneity. Instead two other MR features have been commonly 
used in a large number of experiments to probe the transport sig- 
nature of TSS: the weak antilocalization ( WAL) at low magnetic (B) 
field 19 29 and linear MR (LMR) 30-34 at high B. However, given the 
parallel conducting surface and bulk channels often existing in TIs, 
ambiguities can arise when attributing WAL and LMR to TSS, as 
both phenomena could also arise from bulk states, and have even 
been observed in many non-TI-based spin-orbit-coupled or narrow- 
gap semiconductors 3536 . For example, it is often difficult to exactly 
determine the roles of TSS in earlier WAL studies in TIs with metallic 
bulk conduction, where the reported WAL from only one coherent 
channel likely reflects the strongly mixed contributions from bulk 
and surfaces coupled together 19,24-27 . More recent experiments util- 
izing gate tuning of the surface-bulk coupling and the number of 
phase coherent channels in Bi 2 Se 3 thin films 21,22 have pointed out a 
possibility to extract the transport signatures of the TSS through 
careful analysis. Despite these progresses, most experiments so far 
still dealt with significant bulk conduction, and did not reach or 
clearly demonstrate the so-called "topological transport" regime, 
where bulk is truly insulating and surface dominates the conduction. 
In addition to the ambiguities associated with WAL, the nature of 

LMR 13,30-3436-41 itsdf been debated ( with c l assica p8.39 

and quantum models 40,41 ) since its observation in silver chalcogen- 
ides and other narrow gap semiconductors 36,37 . The recent observa- 
tions of LMR in TIs 3 " 34 have been often interpreted in the framework 
of the quantum LMR 30,32 due to the gapless energy spectrum of TSS, 
although it is not fully clear if the quantum model 40,41 applies without 
a clean separation of bulk and TSS. A better understanding of the 
nature of LMR will benefit its many proposed applications in areas 
such as spintronics or magnetoelectric sensors 32 . 

In this work, we perform a systematic transport study of 
(Bi 0 04 Sb 0 96 ) 2 Te 3 TI thin films (10 nm) grown by molecular beam 
epitaxy (MBE) on SrTi0 3 (STO)(lll) substrates (250 urn, 
Figs. la,b). Using STO (with its very high relative dielectric constant 
£ r at low temperature) as a back gate, we demonstrate a large tun- 
ability of the carrier density (n, by nearly 2 orders of magnitude) and 
Fermi level (Ef) in our film, exhibiting an ambipolar field-effect (FE). 
This allows us to realize a gate-tuned metal-to-insulator transition in 
the bulk of our sample, thus tuning the transport from the bulk- 
dominated regime (where the Ef is in the bulk valence band, BVB) 
to the topological transport regime (where Ef is in the TSS) with 
surface-dominated conduction at low temperature. We also system- 
atically map out the gate and temperature dependent WAL and LMR 
(along with Hall measurements), revealing differences between the 
bulk and TSS transport regimes and providing more insights on the 
nature of such MR phenomena. 

Results 

Bulk metal-insulator transition and ambipolar field effect. The 

temperature dependence of the zero B field longitudinal resistance 
(R xx vs T) measured at different gate voltages V g are presented in 
Fig. lc. The STO substrate, with its very high s r at low T 13,42 , gives a 
strong gate modulation to the sample's n. By increasing V g , Ef can be 
tuned from the BVB to BBG (intercepting TSS, insets of Fig. lc). For 
V g — — 60 V, where Ef is in the BVB, the corresponding R xx decreases 
with decreasing T and saturates at low T, demonstrating a 
characteristic metallic bulk conduction. The temperature below 
which R xx appears to saturate moves to a smaller value at V g = 
— 10 V, suggesting weakened metallic behavior. Further increasing 
V g to — 5 V, R xx shows a clear upturn below —30 K, indicating the 
appearance of an insulting behavior (attributed to freezing-out of 
thermally excited bulk carriers) in the film. Such a bulk insulating 
behavior can be significantly enhanced by further lifting Ef into the 
BBG (V^, = 0 V) and eventually, for V g > 10 V the bulk insulating 



behavior onsets at a T as high as —100 K, with R xx approximately 
saturating with very weak T-dependence for T < 30 K. The bulk 
insulating behavior observed in = 0 V is consistent with the 
ARPES observation of Ef located in the BBG in as-grown films 12 . 
The nearly saturated R xx (terminating the insulating behavior) for 
V g > 10 V indicates a remnant conduction that can be attributed to 
the TSS in the bbg 111618 43 . 44 dominating the charge transport at low 
T (see Fig. SI). Our observations demonstrate a striking transition 
from the metallic to insulating behavior in the bulk of such films, 
driven by Ef (tuned by V g ). This transition can also be regarded as 
that from a "topological metal" to a "topological insulator", and is 
foundational to our study to clarify the relative roles played by the 
bulk and TSS in MR features. 

Fig. Id shows the FE behavior measured at B = 0 T and T = 1.4 K. 
All the curves of R xx vs V g show ambipolar FE. For example, R xx of 
curve "1" is weakly modulated by the gate as V g < — 20 V, but 
increases significantly and reaches a peak of —12 kD when V g is 
increased to —12 V (the charge neutral point (CNP), V CNP , showing 
an on-off ratio of 600%) before decreasing again upon further 
increasing V g . Furthermore, an appreciable hysteresis in R xx vs V g 
depending on the V g sweeping history and direction is observed. This 
hysteresis is common for STO due to its nonlinear dielectric response 
close to ferroelectricity 42 and may also relate to the interface charge 
traps (defects) between STO and TI film. Here, curve "1" (where 
V C np ~ 12 V) represents the first V g sweep from —60 V to 60 V 
after initial cooling down to 1.4 K. The corresponding R xx in curve 
"1" is consistent with the R xx values (shown as boxes with crosses in 
Fig. Id) at each V g extracted from Fig. lc at 1.4 K. However, V CNP is 
shifted to 50 V (curve "2") and 43 V (curve "3") as V g sweeps back- 
ward (from 60 V to —60 V) and forward again, respectively. 
Repeating the V g sweeps from 60 V to -60 V (or -60 V to 60 V), 
the FE curve will stabilize and follow the curve "2" (or "3"), respect- 
ively. All the data presented later are taken after this stabilization and 
in a forward sweeping direction to minimize this hysteresis effect 
(also because of this hysteresis, one should not directly compare the 
Vg's in the following data with those in Fig. lc, instead V g -V C NP ° r 
low- T resistance values are better indicators of the sample state). The 
ambipolar FE in Fig. Id suggests a sign change of dominant charge 
carrier from p-type to n-type as V g crosses Vcot> confirmed by the 
corresponding gate-dependent Hall resistance R xy (exhibiting a sign 
change) and R xx measured at B = — 6 T as shown in Fig. le. We note 
that the charge carriers are holes at V g = 0 V, also consistent with the 
ARPES-measured imposition in as-grown films 12 . We point out that 
the measured V CNP should be understood as an average global CNP 
for the whole sample, as the local CNP (Dirac point, DP, of the 
surface band) position in energy may differ between the top and 
bottom surfaces (as exemplified in Fig. lc inset) and also further have 
spatial fluctuations (leading to electron and hole puddles). Such 
charge inhomogeneity (common in TI 45 as well as in graphene 46,47 ) 
can easily make the minimum sheet conductance (2D conductivity) 
observed near the global CNP (—6 e 2 /h combining both surfaces in 
our sample, Fig. Id) notably higher than the theoretically expected 
minimum conductivity of Dirac fermions 45 47 . 

Gate and temperature dependent magnetoresistance, Hall effect, 
carrier density and mobility. Figure 2 shows the gate-dependent 
AR XX and the corresponding R xy (plotted as functions of B) at 
various temperatures (see also Fig. S2). Here, we define AR XX (B) = 
R XX (B)-R XX (B = 0 T). In Fig. 2(a), all curves of AR^fB) obtained at 
different V g at T = 1.4 K show a gate-dependent cusp at \B\ < 1.5 T, 
a clear signature of the WAL. The amplitude of the cusp can be 
significantly enhanced by varying V g from — 60 V to 60 V to tune 
Ef from BVB to TSS, and reaches a maximum at CNP. Another 
interesting observation is the LMR observed at higher B field in 
AR XX (B), which is also strongly enhanced by gating towards TSS 
and will be discussed in detail later. Meanwhile, the corresponding 
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Figure 1 | (Bi 0 .o4Sb 0 .96)2Te 3 field effect device and its temperature and gate dependent transport, (a) Schematic of MBE-grown (Bi 0 . 04 Sbo.96)2Te3 thin 
films (thickness = 10 nm) on SrTi0 3 (STO, thickness = 250 um, used as back gate); (b) Optical image of a fabricated Hall bar shaped device; 
(c) Temperature-dependent R„ (left axis, with corresponding R„ per square, R n , plotted on the right axis) curves at various different V g . All the curves 
are measured at zero magnetic field (B) and during the cooling process. The upper and lower insets show the schematic band structure with different 
Fermi level [Ej) positions for V g = 130 V and —60 V, respectively; (d) R„ (left axis, with corresponding R n , plotted on the right axis) as a function 
of gate voltage ( V g ) for different V g sweeping histories (all measured at B = 0 T). The arrows label the V g sweeping directions. The symbols (boxes with 
crosses inside) label the resistances extracted from (c) for each V^at 1.4 K; (e) R„and R xy as functions of V s at B = — 6 T and T = 1.4 K. The unshaded 
(shaded) areas mark the V g ranges with p (n) type dominant carriers. 



R xy vs B also shows a strong gate dependence (lower panel of Fig. 2a) 
with two main observations with increasing V g : 1) the slope o{R xy vs 
B initially increases and is followed by a drop as well as a sign change 
while V g crosses the CNP, a direct manifestation of the sign change of 
charge carriers; 2) the corresponding shape of R xy vs B changes from 
linear to non-linear, suggesting a change from one-band to two (or 
multiple) band transport (due to coexisting surface and bulk 
channels of opposite carriers, and/or electron and hole puddles). 
The AR XX and R xy also show significant temperature dependences 
as shown in Figs. 2b-f. As T increases from 1.4 K, the WAL cusp 
gradually weakens and finally disappears at —25 K, where the LMR 
becomes prominent and starts from very low B (<~0.2 T) for most 
of V^'s (Figs. 2c and S3). Further increasing T (>~40 K), AR XX (B) 
becomes parabolic at low B (<2 T), which becomes increasingly 
evident at further elevated T as shown in Figs. 2d-f, with LMR still 
clearly observable at higher B (>2 T). The ambipolar (sign change of 
slope) and nonlinearity (near CNP) behaviors observed in R xy also 
become increasingly evident with increasing T up to —25 K 
(Figs. 2a-c). However, for T > 40 K, the R xy vs B is always linear 
and has no sign change (Figs. 2d-f), indicating one-band behavior 
with p-type carriers. We note that the gate becomes less effective as £ r 
of STO becomes significantly reduced 42,48 at elevated T (see also Fig. 
S4), rendering the ambipolar FE no longer achievable. Our results 
map out a systematic evolution of both R XX (B, V g , T) and R xy (B, V g , 



T), demonstrating the transport properties in such system are highly 
dependent on Ej (modulated by gating) and the temperature. 

Further studies of the field and Hall effects as well as n and mobil- 
ity are presented in Fig. 3. Figs. 3a,b show the temperature depen- 
dences of R xx vs V g (at B = 0 T) and R xy vs V g (at B = -6 T), 
respectively. Consistent with Fig. 2, for Tup to —25 K, we see again 
the ambipolar FE (Fig. 3a) while V CNP increases from 45 V (at 1.4 K) 
to 80 V (at 25 K), related to the decreased e r of STO substrate men- 
tioned above 42 . The corresponding R xy (Fig. 3b) also demonstrates 
the ambipolar behavior up to —25 K, where R xy (initially negative) 
decreases with increasing V g , followed by an upturn and sign change 
as V g crosses CNP. For T > 40 K, both R xx and \R xy \ monotonously 
and weakly increase as V g increases with no indication of ambipolar 
behavior. We extract n and Hall mobility (ii) from R xy and R xx (in 
Fig. 2) using the one-band model at different V g and Tin the regime 
of p-type carriers (mostly from BVB, see below) where a linear R xy vs 
B is observed, and plot the results in Figs. 3c,d. The carrier (holes) 
density (n p ) is -1.8 X 10 14 cm" 2 at V g = -60 V & T = 1.4 K. As V g 
increases and approaches CNP, n p decreases approximately linearly. 
The similar trend is also observed at higher Ts up to 25 K, while the 
slope of rip vs V g decreases significantly with T > 25 K. An effective 
capacitance (C) per unit area of STO can be calculated from the slope 
and C decreases from -290 nF/cm 2 (e r ~ 82000) at 1.4 K down to 
-5.3 nF/cm 2 (s r - 1500) at 200 K (see Fig. S4), consistent with 
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Figure 2 | Magnetoresistance (JR^ = R^BJ-R^B = OT)) and Hall resistance R^ vs B for different V g at various temperatures, (a) T = 1 .4 K; 

(b) T = 15 K; (c) T = 25 K; (d) T = 40 K; (e) T = 70 K; (f) T = 200 K. 



previously observed strongly T-dependent dielectric behavior of 
STO 42 ' 48 . We also note that n p at V g = 0 V (Fig. 3c) decreases with 
increasing T. This is attributed to thermal excitation of «-type car- 
riers, and confirms that the increased V CNP at higher T (in Fig. 3a) is 
mainly due to the decreased STO capacitance. The temperature and 
gate-dependent fi is shown in Fig. 3d. The mobility at 1.4 K increases 
with increasing V^from —50 cm 2 /Vs at V g = —60 V to —140 cm 2 / 
Vs at V g = 30 V. The similar behavior is observed up to 40 K, while 
for T > 70 K, fi becomes — 1 10 cm 2 /Vs and largely V g independent. 
The inset of Fig. 3d shows a summary of fi (in log scale) vs n p for all 
measured T°s, where the data appear to collapse together and can be 

fitted to fi - n 0 e-"r/"° with n„ = 6X 10 15 cm" 2 , ji Q = 133.4 cm 2 /Vs, 
suggesting that fi is mainly controlled by n but not T (up to —200 K) 



in our system. The measured density-dependent mobility may pro- 
vide valuable input for understanding carrier transport and scatter- 
ing mechanism in TIs 49 ~ 52 , important for developing TTbased 
devices. In the case where R xy is nonlinear with B (seen in Fig. 2) 
due to multiple conduction channels and coexisting holes and elec- 
trons, the one-band model will not yield accurate n (in Fig. S5, n 
calculated from such one-band fits starts to deviate from linear In- 
dependence close to CNP). While a multiple-band model 11,17,43-53 can 
in principle be applied to fit the non-linear R xy , we found however, 
such an analysis does not give unique fitting results (yielding signifi- 
cant uncertainties) in our case. Applying one-band model fitting 
(Fig. S5) for our data measured at T = 1.4 K, the lowest carrier 
(electron) density |n| achieved in our sample is —5.5 X 10 12 cm -2 
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Figure 3 | Field effect and gate-dependent carrier density (n p ) and mobility {//) at various temperatures. Temperature dependences of 
(a) longitudinal resistance R„ vs V g measured at B = 0 T and (b) Hall resistance R xy vs V g measured at B = — 6 T; (c) The carrier (hole) density n p and; 
(d) mobility \i extracted using a 1-band model from Fig. 2 at gate voltages where a linear Hall effect is observed. Inset is ji (in log scale) vs n p at 
different temperatures. Dashed line is an exponential fit (ji = f.i 0 e~"f/"° ) to all data points with «o=6x lo' 5 cm -2 , /i 0 = 133.4 cm 2 /Vs. 



(an overestimate for the actual density) at V g = 60 V (where R xy is n- 
type and only slightly non-linear). This value is smaller than the 
estimated maximum electron density (>1.2 X 10 13 cm -2 , see Fig. 
S5 caption for more details) that can be accommodated in the TSS 
before populating BCB, demonstrating that Ef is already located in 
the upper part (above DP) of the TSS at V g = 60 V, T = 1.4 K. This 
also demonstrates that we can reach a regime where the charge 
carriers are mostly from the TSS and band bending near the surface 
is not significant to populate the bulk carriers (otherwise such a low n 
will not be reached). The corresponding n vs V g at low T (Fig. S5) 
demonstrates that we have successfully tuned fif from the BVB, to the 
lower and the upper parts of TSS (through DP), as the V g is increased 
from -60 V to 60 V. 

Gate and temperature dependent WAL and its enhancement near 
CNP with two decoupled coherent surfaces. Now we present the 
gate and temperature effects on the WAL, which is a manifestation of 
quantum coherent transport in the low-_B MR and observed in our 
sample below 15 K. Figure 4a shows the sheet conductance 
correction dGci(B) = Ga(B)-G[j(B = 0T) vs B at various V g 's 
measured at 1.4 K, where Gq = (L/W)/R xx with W and L being 
the width and length of the channel (between voltage probes) 
respectively. While both the bulk and TSS of TI possess strong 
spin-orbit coupling and can give rise to WAL, we have observed 
that WAL is significantly enhanced as Ef is tuned into the BBG to 
suppress the bulk conduction and decouple the top and bottom 
surfaces. The Hikami-Larkin-Nagaoka (HLN) 54 equation (1) has 
been widely applied to analyze AGa(B) due to WAL: 



*«-sH!*|- k (!)] <■> 

where a is a prefactor expected to be -1/2 for a single coherent 
channel, i// is the digamma function, and B^ = hj 4eL^ is a 

characteristic field (with Lj, = y / Dr^, the phase coherence length, 
and phase-coherence time, D the diffusion constant). It is 
found 22 that even with parallel conducting channels (such as bulk 
and surfaces), the total AGo may still be fitted using Eq. (1) in terms 
of an effective a that reflects the inter-channel coupling (with A = 
2 1 a | representing the effective number of coherent channels). Our 
data in Fig. 4a agrees well (up to 2 T) with the HLN fittings (see Fig. 
S6). The extracted values of |a| (a < 0) and at various V^'s and 
temperatures are shown in Figs. 4b-c, respectively. We find that | a\ is 
strongly gate-tunable and exhibits an "ambipolar" behavior (peaks ~ 
1 at CNP) for all T's where WAL is observed, revealing three regimes 
of behavior as the number of coherent channels and degree of inter- 
channel coupling are tuned by the gate: (I) V g < — 10 V, where iy is in 
the BVB, | oc | ~ 0.5 indicates the surface and bulk are fully coupled 
into one coherent conduction channel (A ~ 1); (II) As V g is increased 
(-10 V < V g < +45 V) to lift E f toward the BBG (TSS), |<x| (andA) 
increases, indicating the top and bottom surfaces start to decouple 
from the bulk and each other, toward forming two channels; when V g 
reaches —45 V, where the is close to CNP, | a\ reaches a maximum 
~1 = |— (1/2 + 1/2) | ) with A ~ 2, corresponding to two fully 
separated phase coherent channels (surrounding the bottom and 
top surfaces); (III) Further increasing V g > 45 V to increase the Ef 
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panel (for L v ) are power law fittings (see text). 



above the CNP in TSS and towards BCB, |a| starts to decrease from 1, 
indicating the top and bottom surfaces start to be coupled again. Our 
analysis also suggests that care must be taken when attributing WAL 

to TSS (especially for a 0.5) 19,24-27 in a TI material with both bulk 

and surface conduction. In addition, we find that the phase coherent 
length La also shows interesting gate dependence (Fig. 4c). In regime 
I, in contrast to the weak gate dependence of | a. | , La at low T ( 1 .4 and 
2.8 K) notably increases with increasing V g and peaks at V g = 
— 10 V. Between regimes II and III, L$ reaches a local minimum 
when V g is near CNP, where |a| reaches a peak. 

Figs. 4b-c also demonstrate the temperature dependence of |a| 
and La in different regimes of V g . We see that in regimes II and III, 
|«| is relatively insensitive to T. In regime I, |tx| moderately decreases 
below 0.5 with increasing T when the TSS is coupled to the bulk at V g 
= — 60 V and —10 V, similar to the behavior previously found in 
samples with bulk-dominated conduction and explained as the sup- 
pression of WAL at high T 22 . Fig. 4d shows the temperature depen- 
dences of La at 4 representative V g 's. Previous studies have 
commonly fitted La(T) to a power-law, with the expectation that 
electron-electron scattering would give La proportional to T~ a5 (or 
T~ 0 - 75 ) in a 2D (or 3D) system 44,55,56 . In Fig. 4d, such power-law 
fittings give La proportional to T 038 , T~ ft57 , T~ 0 - 34 , and T~ 0 - 34 for V g 
= -60 V, -10 V, +45 V and +60 V. We note that the fitted 
power-law T-dependences of La in previous experiments range from 
T~ 0 - 24 (Ref. 57 ), T 05 (Refs. 44, 58 ) to T ft7S (Ref. 22) for different TI 
thin films, and ~T~ a3y (Refs. 59, 60 ) or T~°- 5 (Refs. 61, 62 ) for TI 
nanowires, suggesting that the observed power-law can depend on 
detailed material or electronic properties, possibly related to other 
dephasing processes (not just electron-electron scattering), therefore 
a gate-dependence as we see may not be unexpected. In our case, only 



at V g = — 10 V (where La is also the largest, Fig. 4c) we observe La ~ 
T -05 , close to the predicted behavior of electron dephasing due to 
electron-electron scattering in 2D. However, strong deviation from 
this behavior is observed for other V^'s, where La is also shorter (for 
V g — —60 V, La ~ T~ a3S , where only one coupled 2D conduction 
channel exists; for V g = 45 V, where there are 2 decoupled conduc- 
tion channels, and 60 V, both giving La ~ T~ 034 ), suggesting exist- 
ence of additional dephasing processes at these Vg's (such processes 
may be related to other carrier pockets in the valence band for V g = 
-60 V; and electron-hole puddles near CNP for V g = 45 and 60 V). 

Gate and temperature tunable LMR and its enhancement near 
CNP. We now discuss the pronounced LMR observed. In contrast 
to the standard quadratic MR 41 , the observed LMR does not seem to 
saturate in high fields. Figs. 5a-d show the temperature-dependent 
LMR and corresponding R xy at two representative V^'s. The LMR (at 
6 T) in terms of relative MR (AR XX (B)/R XX (B = 0T)) varies from a few 
percent to —30% depending on V g and T. We note the high-B (>2T) 
MR to be slightly sub-linear (super-linear) for T < 25 K (T > 25 K) 
with the 25 K MR being closest to strictly linear (Figs. 5a,candFig. 2). 
Figure 5e shows the slope (k) of LMR (extracted between 3 T and 
6 T) vs V g at different temperatures. When £ / is in BVB (V g <QV),k 
has very little T-dependence and weakly increases with increasing V g , 
whereas k is dramatically enhanced (by as much as 10 times, and 
becomes much more T-dependent) and approaches a maximum 
near CNP as Ef is tuned into BBG (TSS) at low T. Such an 
observation is confirmed by the temperature dependence of k at 
five V g 's plotted in Fig. 5f. For V g = —60 V and — 10 V, where EfS 
are located in the BVB, k has little temperature dependence in the 
measured T range. As £ris tuned into the BBG (TSS), k dramatically 
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increases with the decreasing T and reaches the highest value near 
CNP at T = 1.4 K. Interestingly, as plotted in the inset of Fig. 5f, we 
find that k vs n p (in log-log scale, only including data points with one- 
band tip for holes as those included in Fig. 3c) at different 
temperatures follow a similar trend with k approximately pro- 
portional to n p ~' (except for the data at 200 K), suggesting that the 
carrier density is important to control k, which is significantly 
enhanced as Ef approaches TSS (see Fig. S8). We note that, at a 
fixed n p , k shows little T dependence up to 25 K, while k notably 
decreases with increasing T for T > 25 K (Fig. S8). In the n-type 
regime (close to CNP), the large LMR is accompanied by prominent 
nonlinearity in R xy (Fig. 5d and Fig. 2c). This observation suggests 
that charge inhomogeneity may play an important role in the 
enhanced LMR, as discussed further below. 

Discussion 

Further understanding of gate-tuned WAL. Our gate tunable WAL 
can also be understood 22 in terms of a competition between the phase 
coherence time (z rp , which does not vary strongly with the V g , Fig. S7) 
and the surface-to-bulk scattering time 22 (tsb> which decreases with 
increasing n 63 ), where the effective |a| generally increases with 
increasing z^lx^ 2 as Vg is tuned towards CNP. When the is in 
the BVB, z v (—hundreds of ps, Fig. S7) is much larger than z SB 
(•dps 63 ), resulting in a single phase coherent channel. As Ef is 
tuned into the TSS, Tsb significantly increases due to the reduced n 
and bulk conduction, and ultimately can become larger than t„, 
realizing two-decoupled channels. The weak increase of |oc| at 
higher T in regimes II and III seen in the inset of Fig. 4d may be 
attributed to a decrease in z (f /z SB , which increases the inter-channel 
decoupling, given that both and z v — ul decreases as T increases 
(Fig. S7) while the Tsb should be relatively constant as both R and n 
(Figs, lc & 3c) change little up to —15 K. However, in regime I, |ot| 
moderately decreases below 0.5 at higher T, where t S b is expected to 
be much shorter than z tp in the measured temperature range 22 . Such a 
decrease of |oc| below 0.5 in presence of strong bulk conduction has 
been attributed to WAL getting suppressed when %„ decreases and 
becomes comparable to the spin-orbit scattering time z so at higher 

J22 

Understanding LMR. There have been two common models 
proposed for the LMR, the classical model by Parish-Littlewood 
(PL) 38,39 and the quantum model by Abrikosov 40,41 . According to 
the quantum model 40-41 , a LMR would occur at the quantum limit 
where the applied magnetic field is so large that only one 40,41 or few 36 
Landau levels (LLs) are populated. This condition is more easily 
satisfied in a gapless semiconductor with linear energy- momentum 
dispersion 40,41 . The theory also predicts that AR XX (magnitude of 
LMR) is proportional to 1/n 2 and has no direct dependence on T 
(as long as T remains lower than the energy gap between LLs and the 
Ef). More recently, another model by Wang & Lei, based on the TSS 
and assuming uniform n, relaxes the requirement of extreme 
quantum limit (instead assuming many LLs are filled and smeared 
by disorder) and predicts a LMR with AR XX oc l/« 64 . On the other 
hand, Parish and Littlewood proposed a classical mechanism for the 
LMR, as a consequence of potential and mobility fluctuations in an 
inhomogeneous electronic system, resulting in admixture of Hall 
resistance into R xx on a microscopic level and a LMR 38,39 . The 
classical model predicts that the relative MR = AR XX /R XX (B = OT) 
(thus slope k) should be proportional to a mobility scale fi s = 
max(|/x|,| A/(j), where fi is mobility and is the mobility 
fluctuation, and the cross over field B c (the magnetic field at which 
the MR curve changes from parabolic to linear) is proportional to 1/ 
fis- The previously reported LMR in TIs have been often interpreted 
in terms of the quantum model 40,41 based on the linear dispersion of 
the TSS 30,32 . However, the studies reported so far have not 
systematically measured the dependence of LMR on n, fi and T, 



while such information is important to unambiguously identify 
and distinguish different mechanisms for LMR as discussed above. 
It also remains unclear whether bulk and surface carriers may 
contribute differently to LMR. We find that none of the existing 
models can fully explain our observed LMR. For example, in the p- 
type one-band carrier regime (inset of Fig. 5f), where is in the BVB 
(because the lowest «p (1.7 X 10 13 cnT 2 ) extracted here is higher than 
the estimated maximum hole density (1.2 X 10 13 cm~ 2 ) that can be 
accommodated in the TSS before populating BVB), we are far from 
the extreme quantum limit (with LMR observable at very high LL 
filling factor, eg. >8,000) assumed in Abrikosov's quantum model, 
and the assumptions of linear band dispersion or TSS as invoked by 
Abrikosov 41 or used in Wang-Lei model also do not apply. We have 
plotted the LMR amplitude (by 8R XX = R XX (6T)-R XX (3T), focusing on 
the B > 3T regime where LMR is fully developed) vs n p (Fig. S9) and 
found it cannot be fitted to a single power-law (either l/n p 2 
(Abrikosov) or l/n p (Wang-Lei)) over this density range, but 
rather appears to cross over from a — l/n p behavior for n p > 4 X 
10 13 cm -2 to — l/n p 2 behavior for n p < 4 X 10 13 cm -2 (except for the 
data at 200K). As E f is tuned into the TSS or CNP, the LMR is 
enhanced and shows strong T-dependence (see Figs. 5e,f and 
S9a,b), while concurrently R xy (B) becomes strongly non-linear 
(Fig. 5d, also Fig. 2c) and exhibits a sign-change (carrier type 
inversion), indicating charge inhomogeneity (such as coexisting 
electron and hole puddles) is significant in this ambipolar regime. 
This is at odds with the T-independent LMR predicted in both 
Abrikosov's quantum model and the Wang-Lei model (which also 
assumes uniformly distributed charge carriers), but instead suggests 
that charge inhomegeity (as highlighted in the classical mechanism) 
may play important roles in the LMR. To address the question 
whether the classical model can describe our observed LMR (in 
both BVB and TSS regimes), it is important to examine the 
correlation between k, B c , and fi (Figs. S10-S12). In the BVB 
regime (where p-type carrier fi can be extracted from one-band 
model), we find that k appears to be approximately proportional to 
fi (consistent with PL model prediction if fi s — fi) up to fi — 100 cm 2 / 
Vs, but becomes poorly correlated with fi for higher fi (Fig. S10). Fig. 
SI 1 shows B c as a function of fi, which is qualitatively (B c generally 
lower for larger fi) but not quantitatively consistent with the PL 
model (predicting \IB C to be proportional to ji, if f.i s — fi). 
Furthermore, we note that PL model should predict \IB C to be 
proportional to k (even without direct knowledge of ji s , which 
could depend on Au). We have examined the correlation between 
1/B C and k (Fig. S12) and find that while such a proportionality may 
hold approximately at relatively high T (>40K), it does not hold for 
the full data set (including the 25K data, where LMR is particularly 
pronounced). In any case, our systematic data have revealed the 
following important points: 1) TSS can strongly enhance the LMR; 
2) the charge inhomogeneity also plays important roles in the 
observed LMR, whose behaviors appear to be qualitatively 
captured by the classical model but several aspects are still not 
quantitatively accounted for. A more complete model likely needs 
to take into account both the full band structure (bulk and TSS) and 
inhomogeneity in order to fully explain the observed LMR. Our 
systematic results on the density, mobility, and temperature 
dependences of LMR (Figs. 5 and S8-12) can provide important 
insights for understanding the mechanisms of LMR and key inputs 
to develop a more complete model. 

Methods 

Material synthesis. The high quality (Bi 0 04 Sb 0 96 ) 2 Te 3 (10 nm-thick) thin films 
studied here are grown by MBE on heat-treated 250um-thick insulating STO (111) 
substrates 12 . The schematic of the sample is shown in Fig. la. Previous ARPES 
measurements have demonstrated that the TSS exists in the BBG and the Ej is located 
below the Dirac point, indicating a p-type doping in the as-grown films 12 . 

Device fabrication and transport measurements. The representative device 
structures are denned by standard e-beam lithography (EBL), followed by dry etching 
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Figure 5 | Gate and temperature tunable linear magnetoresistance (LMR). The LMR, AR^/R^fB = OTJ, and the corresponding R xy as functions of 
magnetic field B at (a,b) V g = —60 V and (c,d) V g = 60 V measured at various temperatures (ranging from 1.4 K to 200 K), respectively; (e) The gate 
voltage dependence of the extracted LMR slope (k, extracted between 3 T and 6 T) at different temperatures; (f) The LMR slope (k) vs Tfor different gate 
voltages. Inset shows k vs n p (in log-log scale) at different temperatures for n p values shown in Fig. 3c (one band p-type carriers). Gray band indicates a 
power law with exponent — 1 (k ~ n p ~ l ). The data for T = 200 K can be fitted by ~n p ~ 2 . 



using Ar plasma. The Hall bar electrodes of the devices are fabricated by another EBL 
process followed by e-beam deposition of Cr/Au (5/80 nm). A Cr/Au (5/100 nm) 
film is e-beam deposited on the back of the STO substrate working as a back gate. 
Transport properties are measured by the conventional four-probe lock-in technique 
with an AC driving current of 100 nA at 17.77 Hz. In a typical device (shown in 



Fig lb), the driving current is applied between electrodes "1,4" and the longitudinal 
resistance R xx and Hall resistance R xy are measured between electrodes "5,6" and 
"3,5", respectively. All the measurements are carried out in a variable temperature ( T, 
from 1.4 K to 230 K) cryostat with a magnetic field B (perpendicular to the film) up to 
±6 T. 
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